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00 ! Abstract 

> 
O 

■ It is shown that no set of reasonable approximations leads to a structure 

, 

■ function as a convolution of the PWIA and a FSI contribution. 

Benhar et al have proposed for the structure function of a composite system a form based 
^ , on the Plane Wave Impulse Approximation (PWIA), corrected for Final State Interaction 
51j , (FSI) PI. Although the result has been critically reviewed before B] a recent preprint 



X 
J3 



takes much the same stand and we renew an attempt to clarify the issue. 

Consider the dominant incoherent part of the structure function per particle due to a 
density fluctuation Pqivi) = exp(iqrri) of a selected particle '1' 

S{q,iu) = im/q)<P{q,y) = -(l/7r)Im^$^|pt(ri)G(^ + E° + ze)p,(ri)|$^) (1) 

G{z) = {z — Ha)~^ is the Green's function of the system. q,uj are the momentum and 
energy transfer in an inclusive reaction from which S is extracted and 0(g, y) is the reduced 
response with the energy loss uo replaced hy y = {m/q){uj — g^/2m) Benhar et al 

proposed 

<t)B{q,y) = {q/m)SBiq,uj) = Re y dy'(j)pwiAiq,y - y')RBiq,y') (2a) 
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RB{q,s) = J dye 'y'RB{q,y) = exp J drp2{r,0;r,0)exp(^ixiiq: 

i J drp2{r,0;r,0)xiiq,r;s) 



— *• exp 
Here 



(2b) 
(2c) 



(m/ q)(t)pwiA = SpwiAiq, ^) = J ( ^ ^ — ) ' (3) 

is the response without core recoil in the PWIA in terms of the single particle spectral 
function P. Rb accounts for FSI and is expressed by means of the diagonal two-particle 
density matrix and an eikonal phase Xiiq^ ^) = (^/O') lo ds'V{r — sq). It is totally off-shell 
in coordinate space with finite integration limits instead of — oo, oo, and reflects a particle, 
knocked-out inside the medium which will ultimately not be detected. Finally Eq. ([2cD is 



the Born approximation of (pE]) for weak V. 

We now attempt to derive (^ within a framework, suitable for the discussion of approx- 
imations. It suffices to consider a non-relativistic, infinite system with Hamiltonian, states 
and energies = T + V, $^ and E^^. Assuming local forces 

p\{ri)H{p^, rn)pq{ri) = H{p^, r„) + C{p^, q) (4a) 
>C(Pi, q) = (Pi + qf/2m - pl/2m = q^ /2m + pizq/m, (4b) 

where q = z. With a correspondingly shifted Green's function T{z) = pJ(ri)G(z)pq(ri) (|l]) 
becomes 

S{q,cu) = -(l/7r)Im^$o |r(^ + + 2e)|$° ^ (5a) 
u + E''^-HA{ri;r,) + C{p„q)+te) ' (5b) 

UJ 



+ E\- HA-i{rj) + Tp^ + Vi{ri, r,) + C{p„ q) + lej , (5c) 

with Vi(ri; Vj) = Sj<2V^(r'i — Vj), the residual interaction of '1' with the core. 

In practice one computes the structure function 5*, retaining parts in the argument of F 
and disregarding the remainder 7. The corresponding response is then 

S = Sd + AS^'jd) ^ Sd = Sdfl + Sd,Fsi, (6) 



with some leading contribution 5*^,0 and a Final State Interaction (FSA) part Sd,Fsi- We 
discuss two choices 

-1 
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+ C{p,,q)-Vi{ri;rj)+tej , -fa = E'^ - Ha + Vi{ri;r,) = E'^ - Ha-i - Tp^ (7a) 
+ E^A-HA-l{J) + Tp^+q + ^e) \ = V,{r,;r,), (7b) 



In a) the kinetic energy in Ha resides in the neglected part 7^ and thus suits a large-g 
approximation, where after absorbing q, the recoiling particle '1' has momentum \pi + q| ~ 
Pi.Pj. Its kinetic energy in the appropriate eikonal approximation is £, Eq. (|4ED. For 
the reduced response one has for weak V H 



0a(g,y) = 

X = Xi + X2 = ijn/q) 



1 r°° 
—Re / dse'^y 
111 
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s 



ds'V{r 



P 
s'q) 



P 



sV{r — sq] 



(8a) 
(8b) 



Eqs. (§) give the lowest order terms in a 1/g expansion of in terms of non-diagonal density 
matrices pi,p2- In particular the asymptotic limit (s = sq) 

roc roc r 

Fo{y) = lim (f){q,y) = (47r)"^ / dppn{p) = (47r)~^ / dpp / ds exp[isp]{pi{0,s)/p) (9) 

<}-*°° J\y\ J\y\ J 

The first cumulant corresponding to (|) is 

Mq, y) = T^Re / dse'^y^^^Raiq, s) = Re [ dy'F,{y - y')Ra{q, y') 

ZTT JO J 



2ti 

Ra{q, s) = exp 



P 

^^ P2{r,0;r~sq,0) , ^.^(^^-p.,) 



exp 



pi(0,s) 



pi(0,s) 

In case b) one starts from the exact shifted Green's function except for the residual 
interaction Vi. Insertion of a complete set ^a-i i^^*^ ( |5aD leads to the PWIA 

p + qf 



Sb{q,uj) = SpwiA{q,uj) = J dpp(^p, 
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(11) 



with P{p, E) the single-particle spectral function of the target. 

According to (|^), the above is the dominant part with no FSI left, unless one considers 
the otherwise neglected = Vi{ri]rj). Inclusion of generally, non-diagonal core matrix- 
elements of Vi causes grave complications. It has been suggested to replace Vi{ri; Vj) by an 
optical potential V°^* 



(10a) 
(10b) 



{pMlb\p\n') ^ (p,0|\/i(2;j)|p',0) ^ (p|V°^V) (12) 

However, by definition V°^* replaces only diagonal matrix elements of Vi{ri]rj) and the 
approximation (0) is consequently impermissible. 

With a fast recoiling particle also in b) one is tempted to introduce a Fixed Scatterers 
Approximation (FSA) for Vi. This, however, seems not commensurate with the retention 
of a dynamically active Ha-i- The objection is circumvented if E\ — if^-iO) in (13) 
replaced by an average separation energy (A), or equivalently, if closure is applied to the 
states ^n~^i implicit in P, Eq. (0). Neglecting in addition p{ <^ C one finds 

0f^(g,2/)-0a(g,2/-(A)) (13) 

Except for a small shift, relevant only around the quasi-elastic peak y = 0, closure reduces 
case b) to a). 

We return to the expression (0) of Benhar etal. passing over many intermediate 

heuristic steps in its construction. Its form resembles ([T0|) which without proof has been 
assumed to also hold if the leading asymptotic part Sas{= Fq) SpwiA- The discussion of 
case b) shows this not to be possible, even when assuming (p!2|). 

Next the authors claim (Q) to be the same as used by Silver 0, itself a re-derivation 
of results in but there are obvious differences. We mention the diagonal p2 in (^, 

different from the semi-diagonal in ( |10b|) and also one potential term out of the two in (|8b| ). 



The second one vanishes only for hard-core interactions, using a non-diagonal p2- 

Putting aside the derivation of (@) we next assess the actual difference, comparing 1/q 
expansions of (pa, 4>pwia, which can be shown to have the same asymptotic limit Fo{y), Eq. 
(^). From Eqs. (H) and (^ one finds for the Fourier transforms of the lowest order FSI term 
(f)i{q,y) = {m/q)Fi{y) 

P2ir,0;r- sq,0) 



(14a) 



Palis) = +i f ^^ P2lr',U;r sq,U) _ ^ _ ^^^^ _ 

J p I Jo 

FpwiM = ^nj dpe^P^' l^nWI'd^ - ^n) +^ J dr^^^^^^^^^ J^' ds'V{r - s'q) 

(14b) 



i dr 



p2{r,0;r,0) p2{r,0;r - sq,0) 

' ds V[r — s q) sV[r — sq) 



p Jo p 



The proof leading to the second equation in ( p.4b|) can be found in [|T0| . One notices diagonal 



and semi-diagonal p2 as opposed to exclusively the half-diagonal one in ( |14a| ) . The two 



expressions (|1^ thus differ and so do their first cumulants. Numerical consequences have 
been discussed in a comparison of inclusive cross section of electrons from Fe and Nuclear 
Matter ||^ . In particular for low energy losses there are considerable differences in the relative 
magnitude of the leading and FSI parts, as well as in the total results for ([10|). 



We return to the simplifying assumptions made in the introduction: Should short-range 
repulsion in V produce large integrals in (10c), their contribution may be tempered by a 
re-summation V — > Veff = t (cf. Eq. (^)). Also, the above presentation is not different 
for finite targets, nor is it qualitatively modified when minimal relativistic effects are 
applied. 

We conclude that, starting from a well-defined theory, there seems to be no way to 
actually derive the total response in a form (Q) which modifies a leading PWIA by FSI 
terms. 

I thank S.A. Gurvitz for useful comments and suggestions. 
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